Morphogen Transport in Epithelia 

T. Bollenbach, 1 K. Kruse, 1 P. Pantazis, 2 M. Gonzalez-Gaitan, 2 and F. Jiilicher 1 

1 Max- Planck- Institute for the Physics of Complex Systems, 
Nothnitzer Strasse 38, 01187 Dresden, Germany 
2 Max-Planck-Institute of Molecular Cell Biology and Genetics, 
Pfotenhauer Strasse 108, 01307 Dresden, Germany 
(Dated: February 9, 2008) 

Abstract 

We present a general theoretical framework to discuss mechanisms of morphogen transport and 
gradient formation in a cell layer. Trafficking events on the cellular scale lead to transport on 
larger scales. We discuss in particular the case of transcytosis where morphogens undergo repeated 
rounds of internalization into cells and recycling. Based on a description on the cellular scale, 
we derive effective nonlinear transport equations in one and two dimensions which are valid on 
larger scales. We derive analytic expressions for the concentration dependence of the effective 
diffusion coefficient and the effective degradation rate. We discuss the effects of a directional 
bias on morphogen transport and those of the coupling of the morphogen and receptor kinetics. 
Furthermore, we discuss general properties of cellular transport processes such as the robustness 
of gradients and relate our results to recent experiments on the morphogen Decapentaplegic (Dpp) 
that acts in the fruit fly Drosophila. 
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I. INTRODUCTION 



Morphogens are signaling molecules which are secreted from cells in a restricted source 
region and provide signals to cells located at a distance from this source. They play a key role 
for the determination of cell fates in animal development lj. While the term "morphogen" 
was coined by Turing in his seminal work on pattern formation in reaction-diffusion systems 
P|, the modern paradigm of morphogen action was introduced by Wolpert in 1969 j^. 
According to this paradigm, morphogens spread from the source region into the adjacent 
target tissue where they are partly degraded. The combination of the localized production of 
morphogens, transport, and degradation leads to the formation of a non-equilibrium steady 
state in which the morphogen concentration decreases with increasing distance from the 
morphogen source. This concentration profile is called "morphogen gradient" and contains 
positional information about the distance from the morphogen source. Cells in the tissue 
detect the local morphogen concentration via receptor molecules that are present on their 
surface and respond by expressing a set of target genes in a manner that depends on the 
detected morphogen concentration. In this way, the morphogen gradient can generate a 
pattern of differentiated cells in the target tissue. In the last two decades, the existence of 
morphogen gradients has been supported by considerable experimental evidence. Prominent 
examples of signaling molecules that function as morphogens are Bicoid which acts in the 
embryo of the fruit fly Drosophila , Decapentaplegic (Dpp) which acts in the Drosophila 
wing disk , and Activin which acts in the embryo of the frog Xenopus |a|. 

The mechanisms by which morphogens are transported and gradients are formed are so 
far not well understood. A difficulty in the study of morphogen kinetics is the fact that 
morphogen transport in a tissue is coupled to cellular trafficking processes. It is influenced 
for example by ligand-receptor binding, the endocytosis of ligand-receptor pairs, and the ki- 
netics of receptor numbers. For a long time, it was taken forgranted that morphogens move 
by diffusion in the extracellular space surrounding the cells g In a few cases, there is exper- 
imental evidence for this: the morphogen Activin in Xenopus is an example [1(1] . However, 
for one of the best studied model systems, the morphogen Dpp in the Drosophila wing disk, 
experiments have called diffusive transport into question and suggested an important role of 



cell surface molecules in this process [6 



ll| . Consequently other transport mechanisms than 



L-onsequf 



extracellular diffusion have been suggested p, [Hi, |12| ■ Firstly, Dpp could be transported by 
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transcytosis. Here, transport is achieved via repeated rounds of morphogen binding to cell 
surface receptors, internalization into the cell and subsequent externalization, and release of 
the ligand from the receptor at a different position on the cell surface js, 13|. Secondly, Dpp 
might move by passive diffusion on the cell surface. Here, a certain type of large molecules 
(called Heparan Sulfate Proteoglycans, HSPGs) which are located on the cell surface could 



function as a "carrier" for the morphogens 



Furthermore, morphogen transport could 



occur in cytonemes which are l ong membrane tubes that connect the morphogen source 
cells to cells in the target tissue 12j. Due to the complexity of the problem, a combination 
of theoretical descriptions of morphogen gradient formation and systematic experiments is 
needed to identify the dominant morphogen transport mechanism jllll . Il^ . 

During animal development, the precision of the positions of differentiating cells in the 
tissue and the times at which cells differentiate is typically high 2 5 Q|- This indicates 
that robust mechanisms that are insensitive to changes of environmental conditions and to 
intrinsic fluctuations have evolved to control cell differentiation. Clearly, such robustness 
could be achieved if morphogen gradients are themselves robust. Evidence for the robustness 
of morphogen gradients was found in recent experiments flfiL . This robustness must 
originate in the mechanisms by which morphogens are transported and degraded. The 
robustness and precision of morphogen gradients 3, 18, [3, 



20, 



2l( or a possible role of 



22j can only be understood using a combination 



anomalous diffusion in morphogen transport 
of theoretical and experimental efforts. 

In this article, we provide a detailed description of morphogen gradient formation by 
different mechanisms and provide full derivations of the morphogen transport equations. 
Some of our results have been recently published in j^j]. Our description captures several 
processes that are supported by experimental data for the morphogen Dpp in the Drosophila 
wing disk. These experimental findings will be briefly summarized in section |H] While we 
focus on the wing disk of the fruit fly here, these processes are very likely to play an important 
role for morphogen gradient formation in other animals. Starting from a description of 
cellular trafficking processes, we derive in section UTTl effective transport equations on larger 
scales. We investigate key properties of gradient formation by these mechanisms in section 
IIVI and study the effects of a directional bias for transport resulting from cellular polarity. 
Finally, we extend our approach to higher dimensions in section IS and discuss morphogen 
transport in two dimensional epithelia. 
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II. MORPHOGEN GRADIENT FORMATION AND CELLULAR TRAFFICKING 
PROCESSES 



le so-called wing primordium. In the 



13. 



14j . Dpp is produced and secreted 



The larva of the fruit fly Drosophila contains precursors of the organs of the adult animal. 
The precursor of the fly wing is a flat pouch that consists of two cell layers that are connected 
at the edges and is called wing disk (see Fig. 1 in |l3j]). The thicker one of these cell layers 
is formed by columnar epithelial cells and includes t 
following, we consider this two dimensional cell layer 
in a specific source region which is a narrow stripe with a width of about 7 cell diameters 
that is located at the center of this layer. Cells outside of this source do not produce Dpp 
but posses receptors located at their cell surface to detect its presence. Dpp spreads from 
the source region into the adjacent target tissue on both sides of the source region. In 
the whole tissue, Dpp molecules are degraded. As a consequence of the localized source 
and degradation, a graded morphogen profile is built up. This formation of the morphogen 
gradient can be directly observed in experiments by using a Dpp that is labeled with green 
fluorescent protein (GFP-Dpp) 0,Q|. In steady state, the Dpp gradient extends over 50/mi 
into the target tissue. This corresponds to about 20 cell diameters. 

Several cellular processes are relevant during the formation of this morphogen gradient. 
Morphogens are ligands which bind to specific receptor molecules. Ligand-receptor pairs are 
internalized into the cell via endocytosis. Internalized ligands are either degraded or they 
can be recycled to the cell surface via exocytosis. Finally, receptor-ligand pairs can diffuse 
on the cell surface and free ligands can diffuse in the extracellular space surrounding 
the cells. 

Furthermore, cells in the wing disk produce, degrade, externalize, and internalize re- 
ceptor molecules. In general, the production rate of these receptors is affected by the local 
morphogen concentration. For example, a high Dpp concentration leads to a reduced produc- 
tion rate of the Dpp receptor Thick- veins (Tkv) in the wing disk |2( . This behavior is called 
"receptor down-regulation". The opposite behavior can also occur: high concentrations of 
Hedgehog (Hh), another mor phog en acting in the wing disk, lead to an increased production 
of its receptor Patched (Ptc) 24j. This phenomenon is called "receptor up-regulation" . 

Recent experiments have revealed the interplay of Dpp gradient formation and cellular 
trafficking. Firstly, endocytosis has been blocked in the whole wing disk except for the 
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Dpp source region or, alternatively, in smaller patches of cells (clones) in the tissue jy, 
Il3| |. This has been achieved using mutant flies in which endocytosis can be blocked at an 
elevated temperature of 34°C due to the temperature-sensitive mutation shibire. Five hours 
upon blocking endocytosis in the whole target tissue, GFP-Dpp fluorescence was almost 
undetectable in the target tissue while its gradient extended visibly over more than 20 cell 
diameters into this tissue before the endocytic block. When endocytosis was only blocked 
in patches of cells near the Dpp source region, a pronounced transient depletion of the 
GFP-Dpp concentration behind these clones was observed. Both experiments indicate a 
role of endocytosis in Dpp transport. This suggests that instead of simply diffusing in the 
extracellular space, Dpp is transported via the cell interior in repeated rounds of endocytosis 
and exocytosis. This transport mechanism is called transcytosis [6(. Indeed, a theoretical 
description in which morphogen transport is solely based on extracellular diffusion while 
interactions of the morphogen with its receptors are taken into account was found to be 
inconsistent with the experimental observations Q|. 

In a different set of experiments, the role of large cell surface molecules (HSPGs) in 
morphogen transport was investigated ^J. A depletion of extracellular Dpp was observed 
behind mutant clones of cells which lack HSPGs. This indicates a role of HSPGs in Dpp 
transport. Indeed, it has been suggested that HSPGs facilitate the diffusion of morphogens 
on the cell surface 



III. MORPHOGEN TRANSPORT IN ONE DIMENSION 



We introduce a discrete description of morphogen transport by transcytosis and passive 
extracellular diffusion and derive transport equations for these processes. The theoretical 
description developed here is generally applicable to biological systems in which molecules 
are transported by extracellular diffusion and transcytosis. To stress this generality, we will 
mostly refer to the transported molecules as "ligands" instead of "morphogens" . 



A. Ligand kinetics in a chain of cells 

In one space dimension, we describe the ligand kinetics in a chain of cells, see Fig. We 
denote the distance between the centers of two neighboring cells by a and the width of the 
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gap between two cells by b. In this chain, n indexes the cells, see Fig. ^ The number of free 
extracellular ligands between cells n and n + 1 is denoted L n . The numbers of intracellular 
free and ligand-bound receptors are denoted R$ and S$ , respectively. R$ and R$ denote 
the concentrations of free receptors on the left and right cell surfaces, respectively. S^' and 
S^ r ) denote the ligand-bound receptors accordingly. The kinetics of the ligand and receptor 
numbers are given by 



d L 
It 1 

dt n 

d Q ( r 

dt n 
dt n 
dt n 
dt n 



kos(S^ + S'n+i) ~~ k on {R^ + R^ +1 )L n H ^(L n+1 + L n _i — 2L n ) — e^egL 

cr 

f s V n i U a oO) _ U E>(r) T _ f. D(r) , 

—k gS^ ^ + k on R n L n — b- m tS^ n ^ H 2~^n 

fextRfJ ~l" f'mt{RrJ -^n ) fdegRn^ 



(1) 



Here, the binding and un-binding of ligands to and from receptors is characterized by rates 
fcon and k ff. The internalization and externalization of receptor-ligand complexes is captured 
by the rates 6j nt and 6 ext . Free ligands can hop directly from one gap between the cells to the 
adjacent ones at a rate 2D /a 2 . This describes their free diffusion in the extracellular space 
around the cells with diffusion coefficient D . The degradation of ligands in the extracellular 
space occurs with rate ed eg and that of ligands bound to receptors inside the cell with rate 
bdeg- Furthermore, free receptors are internalized and externalized with rates fi nt and f ext 
respectively. Internalized free receptors are degraded with rate /d eg . In addition, each cell 
produces receptors with a rate f syn . 

Newly produced receptors appear on the cell surface [25]. The rate of receptor synthesis 
/ syn in Eq. © depends on R® + R$ and S® + S®: 



/syn f t 





syn 



Ru 



(2) 



where / s ° yn is a basal rate of receptor synthesis and -R max is the saturation value of the surface 
receptor concentration at which the production of new receptors stops. The dimensionless 
parameter ip couples receptor synthesis to the concentration of morphogens. The two cases 



of receptor up- and down-regulation 



25| are captured by ip < 1 and ip > 1, respectively. 



The description (0) is valid in the bulk of the system. We still have to specify the kinetics 
at the boundaries of the chain. To describe the effects of a ligand source located at n = 
in Eq. (JTJ, we modify the equation for the free ligand: 

— L = k gS® — k on R^pL + —^(L\ — L ) — e^gLo + u, (3) 

where v is the rate at which ligands from the source enter the system. At the position n = N 
where the lattice ends, we impose 

— Ln = k gS$ — k on Rpf Ln + — ^(Ljv-i — Ln) — edeg-^iv, (4) 

which describes a zero flux boundary condition at this edge of the lattice. This boundary 
condition can be motivated by the geometry of the Drosophila wing disk Q]. However, if 
N is sufficiently large the ligand number at N is small and the choice of boundary condition 
has only a small influence on the ligand profile. A sequence of morphogen and receptor 
profiles at different times that were obtained by numerical solution of Eq. (0) are presented 
in Fig. 121 together with the steady state profiles that represent the morphogen gradient. 



B. Effective transport equations on larger scales 

We derive effective continuum transport equations for ligand transport starting from Eq. 
(JTJ) - We introduce the concentrations l(t, x) = L n (t)/a, r^ l \t,x) = R$(t)/a, r^ T '(t,x) = 
R$(t)/a, ri(t,x) = R®{t)/a, s®(t,x) = S®{t)/a, s^(t,x) = S®(t)/a, and Si(t,x) = 
S$(t)/a where x = na. We derive continuum equations for the kinetics of these densities 
starting from Eq. (JIJ) by locally expanding the densities in a power series with respect to x, 
for example L n+1 /a = l(x + a) = l(x) +ad x l(x) +a 2 9^/(x)/2. It is further useful to introduce 
the new variables r±(t,x) = r^'(t,x) ± r^ r \t,x) and s±(t,x) = s®(t,x) ± s^ T '(t, x) so that 
r + and s + measure the total free and ligand bound surface receptor concentrations per cell 
and r_ and s_ the polarization of these concentrations on the cell surface, respectively. 

In situations where the length £o over which the steady state gradient decays is large 
compared to the cell diameter a, a separation of time scales occurs in the system which 
makes the adiabatic elimination of rapid variables possible. Indeed, if r a is the relaxation 
time of the kinetics within one cell, the slow relaxation of the gradient occurs on a time 
scale t^ d = r a (££>/a) 2 ^> r a . We thus use the approximation that all local kinetics relaxes 
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instantaneously. At each position x, this yields the relations 




h nlr — r f) 1 -i- ("-"jQcxt Q 



s 




(5) 



Using these expressions, we can adiabatically eliminate five of the seven variables 
I, Si, s+, s_, Tj, r + , r_ and obtain only two coupled equations for the remaining slow vari- 
ables which are the total ligand density X(x, t) = l(x,t) + Si(x,t) + s + (x,t) and the total 
receptor density p(x, t) = r^x, t) + r + (x, t) + Si(x, t) + s + (x, t): 



The other densities I, s i} s + , s_, r i} r + , r_ can be calculated from A, p, and their first spatial 
derivatives via Eq. (jHJ). The derivation of Eqs. ©, ©> and (J7J) is discussed in appendix|X| In 
addition to equations © and ((7j), this procedure provides us with explicit expressions for the 
effective diffusion coefficient D\, the effective degradation rate k\, the receptor degradation 
rate k p as well as the transport coefficient D p which describes ligand transport induced by 
gradients of the receptor concentration, see Eq. (|A6|) . In Fig. El A,B, the coefficients D\, 
D p and k\ are displayed as a function of the ligand concentration A for a typical choice of 
parameters. The effective receptor production rate z/ syn (A,p) = f syn /a is a function of the 
ligand and receptor concentrations. The functional form of u syn that corresponds to the 
receptor production rate (J2J) in the discrete description is 



with r max = -Rmax/cJ- Ligand transport is associated with the current j = — (D\(\, p)d x \ + 
D p (\, p)d x p). The current j at x = is related to the secretion rate v in Eq. (j3J). In one 
dimension jo = v. 

In Fig. |21 we show time-dependent receptor and ligand profiles which are solutions to 
Eqs. (0) and (JJJ). These solutions are in agreement with the corresponding solutions to the 



d t \ = d x (D x (X,p)d x X + D p (X,p)d x p)-h(X,p)X, 

dtp = ^syn(A,p) - k p (X,p)p. 



(6) 
(7) 
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discrete description (J!} which validates the adiabatic approximation made in the derivation 
of Eqs. © and ©. 

C. Effects of a directional bias 

If cells possess a polarity, transcytosis can have a bias and lead to directed transport. 
Here, we include a directional bias in our description of ligand transport. Experiments 
indicate that the transport of the morphogen Dpp in the Drosophila wing disk is non- 
directional on macroscopic length scales However, epithelia with cell polarity could in 
principle exhibit directional transcytosis. 

We therefore generalize the discrete description (0) by allowing receptor-bound ligand 
molecules to be preferentially externalized on either the left or the right cell surface. To this 
end, we introduce a dimensionless parameter j3 which measures this bias: for f3 = —1/2, all 
receptor-bound ligands are externalized on the left surface, for (3 = 1/2 on the right surface, 
and for (3 = 0, we recover the unbiased description, see appendix |B] for details. 

In presence of such a bias, the transport equations © and (J7|) generalize to 

d t X = d x (D x (X }P )d x X + D p (X,p)d x p-V (X,p)X)-k x (X,p)X, 

dtp = u sya (X,p) -k p (X,p)p, (9) 

where Va(A,p) is a concentration-dependent effective drift velocity, see Eq. (jB3|) . The other 
coefficients in Eq. Q remain the same as in the case without directional bias. In Fig. 0] A, 
we show V/3 as a function of A for (3 = 0.1, i.e. a bias that leads to preferential transport to 
the right. Typical ligand and receptor profiles that are generated in presence of a bias are 
shown in Fig. 0]B,C. 

D. Constant surface receptor concentration 

We now discuss the simple case where the total surface receptor number R is constant 
everywhere. This approximation is useful because it still captures most important features 
of morphogen transport by transcytosis j^J. 

In this case the receptor kinetics in Eq. becomes obsolete. The equations describing 



9 



ligand transport in this simple case are 

— L n = k ff(S^ + S'n+l) ~~ kon(R — S$ — Sn+i)L n + —^-(L n+ i + L n -\ — 2L n ) — edegL n 

:Si ■* = — k Q ftSi ^ + k on ( — Si ^)L n — b- mt Si ■* + —b ex tSi^ 



^ n '"on"n ' '-on y 2 n ) » "mi~n 1 2 



^Sg> = -fc^Sg) + 6 int (5« + Sf>) - 6 deg 5«. (10) 

Boundary conditions analogous to Eqs. and (J3j) are imposed at n = and n = N. 

If the surface receptor concentration is constant, the continuum limit after adiabatic 
elimination of fast variables is described by 

9 t X = d x (D(X)d x X) - k(X)X , (11) 

where the effective diffusion coefficient D(X) in the absence of extracellular diffusion (for 
D = 0) and the effective degradation rate k(X) are given by 

a 2 b ext b int k oS ak on rC - (A) 



D(X) 



4:A(X)(2ak on rk o{i (b ext + b int ) + 6 in tC_(A)) 



= ( b deg b int e deg fc off \ 

1 ; ~ ak on X \2b ext (b ext + b int ) C_(A) ) ' 1 1 

In these expressions, r = R/a and 

r i V 2 

^(^) = [-4&cxt(&cxt + b- int )a 2 k 2 on rX + (b hlt ak on r + b cxt B + (X)) 2 

B±(X) = k oS + ak on (X ± r) 

C±(A) = b int ak on r =p A(A) ± 5 cxt £? ± (A). 

These coefficients -D(A) and fc(A) are shown as a function of A in Fig. EjC Their non-linear 
dependence on the total ligand concentration A is very similar to that of the coefficients 
D x (X,p) and k\(X, p) in Eq. ®, see Fig. El A. 



We can describe transport by diffusion of receptor-bound ligands in the cell membrane 
by the same methods. We discuss this mechanism in appendix [D] where we consider the 
case where endocytosis and recycling are unimportant. Furthermore, we have so far dis- 
cussed the case where extracellular diffusion is weak (D Q small). The opposite case in which 
extracellular diffusion dominates ligand transport is discussed in appendix IEI 
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IV. PROPERTIES OF TRANSPORT BY TRANSCYTOSIS 



A. Nonlinear diffusion and degradation 

Several key features of transcytosis follow directly from the general shape of the transport 
equations (0) and (J7|) and from the dependence of the transport and degradation coefficients 
on the ligand and receptor concentrations. Firstly, the presence of the term D p (X, p)d x p 
in Eq. © shows that gradients of the receptor concentration induce a ligand current that 
is directed towards higher receptor concentrations since D p < 0. This contribution to 
the current comes up, because the ligand affinity for a region increases with the receptor 
concentration in that region. For small ligand concentrations D p ~ A which ensures that the 
corresponding current vanishes. Furthermore, D\ and k\ in Eq. (JBj) approach finite values 
for small A. In this limit, Eq. © consequently becomes a linear diffusion equation with 
degradation. This implies that on large length scales and for small ligand concentrations, 
transcytosis is indistinguishable from passive diffusion. 

In the opposite limit of large ligand concentrations A, D\ and k\ exhibit the asymptotic 
behavior 

Da- A) + Cl (p)A- 2 

and k\ ~ ed cg + C2(p)A _1 . Here D is the extracellular diffusion coefficient which is ap- 
proached in the limit of large A, and we have defined ci(p) = afeext/int&oftp/4/c n(fee X t + feint) 
and c 2 (p) = fedeg feint p/(feext + feint)- The transport coefficient 

D p ~-D - Cl (p)\- 2 

in this limit. Interestingly, both D\ and \D p \ approach the value D Q . This implies that 
transport is dominated by extracellular diffusion for large ligand concentrations A. This 
behavior results from the fact that most receptors are occupied and can consequently not 
participate in ligand transport by transcytosis. A maximum of D\ can occur for intermediate 
values of A as long as Dq is smaller than a critical value, see Fig. |3] A. Similarly, there can 
be a minimum of D p as a function of A, see Fig. EJB. The observation that ligand transport 
is most efficient at a specific ligand concentration A is an interesting consequence of the 
nonlinearities of the transport process. 

In the special case D Q = 0, D p as well as D\ vanish in the absence of receptors, i.e. 
for p = 0, or if either binding or un-binding of ligands from the receptor, internalization 
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or externalization of occupied or free receptors is suppressed, i.e. if either one of the rates 
k on , k s, hnt, b ext , fmt, or f ext vanishes. This reflects that in the absence of extracellular 
diffusion, transport is generated by repeated internalization and externalization of ligand- 
bound receptors as well as ligands binding to and un-binding from surface receptors. In the 
limit of fast internalization or fast un-binding, the ligands are confined to the cell interior 
or the extracellular space, respectively, and transport is consequently hampered. Indeed for 
D = 0, D\ — > if 6; n t — > oo or if k Q s — > oo. We discuss other limits of our description in 
appendix IU1 

Similarly, if transcytosis has a directional bias, the effective drift velocity Vp vanishes if 
either p, k on , k o{i , b int , b ext , f int , or / ext is zero and also for b int -> oo or k oH -> oo. Moreover, 
the drift velocity Vp is independent of D . For small A, Vp adopts a finite value. \Vp\ 
is a monotonically decreasing function of A and, in the limit of large A, Vp vanishes as 
Vp ^ c 3 (p)A~ 2 , with c 3 (p) = pb cxt f- mt k ofi p/[(b ext + b int )k on ], see Fig. H A. The fact that Vp 
vanishes asymptotically for large A again reflects that transport is mediated by receptors 
which are only present in limited numbers. 

B. Steady state concentration profile 

We now calculate the steady state ligand profile formed in the half space x > in the 
presence of a source which is located at x < 0. In the steady state, Eq. (J2j) with the 
condition d t p = yields a relation p(x) = p s (X(x)) between the total receptor concentration 
p and the total ligand concentration A at position x. This combined with Eq. (JHJ) leads to 
the steady state equation for the ligand profile 

d x (D s (X)d x X)-k s (X)X = 0, (13) 

with the effective diffusion coefficient in the steady state D S (X) = D\(X, p s (X)) + 
D p (X, p s (X))dp s (X) / dX and the effective degradation rate A; S (A) = k\(X,p s (X)). The steady 
state relation p s (X) is a monotonic function of A and converges to finite values for A — > oo 
and po for A — > 0. This reflects that each cell only contains a limited number of receptors 
and is never completely devoid of receptors. It implies that dp s /dX = for large A, so that 
in this limit D a (X) ~ L>a(A,Poo)- 

The steady state ligand profile X(x) described by Eq. (fT3*j) can be determined exactly. 
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We rewrite Eq. (JT3j) as 



d x j s = -ks(X)X 

d r X - ux) 



x 



[ X(X) dX'D s (X')/UX') , (15) 
Jx(o) 



A (A)' 

where the steady state current j s is a function of X(x) only. This implies 

djs(X) k s (X)XD s (X) 

dX j e (\) 
dx{X) D S (X) 

dX " Js (A)' 1 1 

where x(X) is the inverse function of the steady state ligand profile A (a;). Using Eq. (JHj), 
we find the steady state solution 

r x(x) 

/A(0) 

where the steady state current is 

j 8 (X) = (2 £dX' k s {X')D 8 (X')X'\ ' . (16) 

In the steady state, the total ligand concentration decreases monotonically with increasing 
distance to the source. For small A, the ligand profile decays as A ~ exp(— x/£) with 
£ = \jD s (0)/k s (0). For large ligand concentrations A ^> At and in the absence of free 
diffusion (i.e. for D = 0), the current j B behaves asymptotically as 

j s 2 (A) ~ j's(A T ) 2 + 2e dcg c 1 (p 00 ) ln(A/A T ) + 2c 1 (p 0O )c 2 (p oo )(l/A T - 1/A). 

Here, A^ denotes a crossover value beyond which the asymptotic behavior becomes valid. 
Therefore, the current diverges logarithmically as j 2 ~ 2ci(p 00 )ed eg In A for large A and 

e dcg > 0. 

This behavior of j s has interesting implications for the steady state ligand concentration: 
X(x) is characterized by a singularity which occurs at a position x* < that moves towards 
x = as A(0) becomes large. In the vicinity of x*, X(x) behaves as 

A~ (x-x*)-\-]n(x-x*))-^ 2 . (17) 

Note, that the case ed eg = has to be discussed separately. In this case, the current reaches 
for large A a finite maximal value j max and the steady state profile diverges as 

A ~ ci(p 00 )/(x - a;*)j m ax- (18) 
13 



If extracellular diffusion is present, i.e. if Dq > 0, D S (X) in Eq. ()13|) changes its asymptotic 
behavior to D s ~ c^p^/X 2 + D . For large A > Xd with Xd — (ci(poo) / A)) 1//2 , Eq. 

becomes linear and the steady state solution decays exponentially on a length scale 
£d = \J Dq/ e^cg • The nonlinear behavior described by Eqs. (|17|) and (J 18)) is thus valid for 
Xx < A < A/j. 

C. Robustness of morphogen gradients 

To study the robustness of morphogen gradients with respect to changes of the morphogen 
secretion rate, we consider the response of the steady state gradient to changes of jo- We 
define the following dimensionless measure of robustness: 

K{]q,X) = a{j d jo x{\))-\ (19) 

where x(X) is the position at which the steady state ligand profile attains the concentration 
A. Here, a robustness of 7Z{X) = 1 implies that under a 100% increase of jo the position at 
which the ligand profile attains the fixed value A is displaced by about one cell diameter a, 
see Fig. 03 Thus for 7Z(X) > 1, the shift of the position x where the ligand concentration 
has the value A cannot be detected by the cells in the target tissue even under significant 
changes of j . 

In absence of extracellular diffusion, the singular behavior ()17J) of the steady state profile 
near x = x* has remarkable consequences for the robustness of gradient formation. Using 
the robustness 7Z defined in Eq. ffT9")) . we see from Eq. ([TB]) that 1Z is independent of A. For 
steady state equations of the general form (fTT?J) . the robustness can be calculated as 

TZ = a(j d jo x)~ l = ad Xo jo/D s (X ) = ak s (X )X /jo (20) 

where A = X(x = 0) and Eqs. ([To]) and (fTB)) have been used. The robustness is thus 
completely determined by the ratio of the effective degradation rate and the ligand current 
at x = and does not depend on A, i.e. it is the same at all positions x of the concentration 
gradient X(x). High degradation rates and small currents lead to a robust gradient. Using 
the asymptotic behavior of the steady state profile for D = 0, we find that the robustness 
increases rapidly for large currents jo as 1Z ~ Jq 1 ^^^ with j 2 = l/2ci(poo)edeg- For small 
j , 1Z ~ a/£ becomes constant. In Fig. El A, B, we illustrate the behavior of the robustness 
of steady state gradients for a small and a large value of j . 
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The situation is different if free diffusion in the extracellular space is present. As discussed 
in the previous section, the singularity in the steady state solution disappears for D > 0. 
As a result of this, the robustness approaches a finite value 7?. max = a/£d as j — > oo. In 
Fig. |5]D, 7Z(jo) is shown for different values of £d/a and Fig. shows an example for the 
effect of the presence of extracellular diffusion on the robustness of the gradient. 

In summary, we find that morphogen gradients can be extremely robust to changes in 
the morphogen secretion rate of the source cells if transport is dominated by transcytosis. 
The presence of extracellular diffusion reduces this robustness that is completely lost when 
extracellular diffusion is the dominant transport mechanism. 



V. MORPHOGEN TRANSPORT IN TWO DIMENSIONS 



A. Ligand kinetics on the cellular scale 

The theoretical framework introduced in the previous sections can be extended to higher 
dimensions. Considering a two dimensional geometry, we represent the cells in the tissue on 
the sites of a discrete lattice. The tissue geometry of the wing disk can be captured by an 
irregular tiling of the plane, see Fig.|U]A,B. For simplicity, however, we use a triangular lattice 
with N hexagonal cells in our description, see Fig. 13 D. Each hexagonal cell n = 1, . . . , N 
has J edges (n, j) with j — 1, . . . , J (J = 6 for a triangular lattice) along which it faces a 
uniquely defined neighboring cell n' at its edge (n',j f ), see Fig.|H]C. The space between these 
two facing edges of the neighboring cells is denoted by the symbol (n, j) with the property 
(n,j) = (n',f). 

To keep our notation simple, we first discuss the case of constant surface receptor con- 
centration, see section IlII Dl Assuming that there is no directional bias of ligand transport, 
the externalization of the ligand-receptor complexes can occur on all surfaces of the cell 
with equal probability. The equations for the ligand kinetics corresponding to Eq. (fTU|) on 
a lattice are 

d . , ,„ „ s Jk nn /2R 



L(n,j) — k s(S n j + S n ' t j>) 2 ( j S n j S n ij^j L( n> j) edcg-^(nj) 

S n ,j = — (feint + k fi)S n ,j H ^— (-J — Snj^j 



dt 

d -S$ = -b cxt S^ + b int J2 - b deg S£\ (21) 



dt n 
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where S$ is the internal bound ligand concentration in cell n, S n j the surface bound ligand 
concentration on edge j of cell n, and L/ n j\ the free ligand concentration in the extracellular 
space (n,j) which is located between the two adjacent cells n and n' and the edges (n, j) 
and (n',j'), see Fig. El C,D. In Eq. (|2ip. we have for simplicity neglected extracellular 
diffusion which would couple the concentration L/ n j\ to the one on neighboring sites. At the 
boundaries of the lattice, the equation for J^L( n j) in (121(1 takes into account ligand influx 
analogously to the one dimensional situation, see Eqs. (JHJ) and (jlj). 

B. Transport equations on larger scales 

The effective behavior of ligand transport as described by Eq. (|21|) exhibits anisotropy on 
large scales due to the anisotropic lattice structure. We consider for simplicity an isotropic 
continuum limit. This simplification is motivated by the irregular arrangement of cells in a 
tissue which does not exhibit lattice anisotropies. 

In the isotropic case, the continuum limit describing transport by transcytosis in two 
dimensions is of the general form 

d t \ = V • 0D 2d (A)VA) - £; 2d (A)A (22) 

The effective coefficients D 2d (X) and fc 2d (A) in Eq. ()22|) are in general different from those in 
the one dimensional case p2|l . In order to determine values for the coefficients _D 2d (A) and 
A; 2d (A) we first consider concentration profiles which vary only along one symmetry axis of 
the lattice given by the x-axis in Fig. El D. In this situation the problem can be represented 
on a one dimensional lattice similar to the one dimensional chain discussed above. We thus 
determined _D 2d (A) and k 2d (X) along this lattice axis using our one dimensional approach. 
Lattice symmetry implies that these coefficients apply to three different lattice axes. In our 
isotropic simplification, we assume that they apply to all directions. From this argument, 
we find the same effective degradation rate k 2d (X) = k(X) as in one dimension, see Eq. (|12jl. 
The effective diffusion coefficient changes by a factor of 2/3: -D 2d (A) = (2/3)D(A). 

By the same considerations, the more general description with receptor kinetics intro- 
duced in one dimension in Eqs. © and (J7J) generalizes in two dimensions to 

d t X = V • (D 2d (A,p)VA + D 2d (X,p)Vp) - k 2d (X,p)X, 
dtp = v™(X,p)-k?(X,p)p. 
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Here the coefficients L>f d (A,p), D 2d (X,p), /c^ d (A,p), and k 2d (\,p) are in general modi- 
fied due to the lattice geometry. For a triangular lattice D 2d (X,p) = (2/3)D\(\, p) and 
D 2 p d (X,p) = (2/3)D p (\,p). The degradation rates k 2 x d and k 2d are identical to those in the 
one dimensional case. Finally, the rate of receptor synthesis u 2d n = u Byn /a is the same as in 
one dimension, see Eq. (JSJ), but measured in different units. At the boundary line at x = 0, 
the ligand source is described by a current jo = v/a across this boundary line. 

We have compared solutions of the effective continuum equation in two dimensions ([22)1 
to those of the discrete description on a triangular lattice fl2*T]) . Here, the ligand source 
located at x < extends along the ^/-direction. This setup is translation invariant along this 
direction if no inhomogeneities are present in the tissue. In this situation, the solutions of 
the discrete and continuum descriptions are in good agreement. In order to test the validity 
of the continuum description for a case that is not translation invariant along the ^/-direction, 
we compared the solutions of Eq. ()22|) for a geometry in which a rectangular region which the 
ligand cannot enter is present in the tissue to solutions of Eq. (}2T| where such a rectangular 
region is approximated, see Fig. In the discrete description, we imposed this constraint 
by setting b- mt = in this region. In the continuum description ()22j) . this corresponds to 
D 2d = within this region which we realized by imposing a zero flux boundary condition 
(j = 0) at its border. As in one dimension, the continuum description is appropriate as long 
as the degradation rates and the hopping rate D /a 2 which describes extracellular diffusion 
remain small. 



VI. DISCUSSION 



In this article, we have first presented a description of morphogen transport in which 
cells are discrete entities. This description is based on key processes like the diffusion 
of morphogens in the extracellular space, binding and un-binding of the morphogens to 
and from receptor molecules that are located on the cell surfaces, internalization of these 
receptor-ligand complexes into the cell and their subsequent recycling, as well as degradation 
of external and internalized ligands, see Fig. ^ Moreover, the production and intracellular 
trafficking of free receptor molecules by the cells is included. We have derived effective 
nonlinear transport equations (jBJ) and ((Zj) for the total morphogen concentration and the 
total receptor concentration which describe transport by transcytosis on larger length scales. 
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The effective diffusion coefficient and the effective degradation rate in these equations are 
concentration-dependent. If transcytosis has a directional bias, an additional drift term 
appears in the transport equations. 

Other mechanisms of ligand transport can be effectively described by equations (jUJ) and 
(J7J). The effective transport coefficients can be derived from a detailed description of any 
particular transport mechanism. As an example for this, we have discussed a model of mor- 
phogen transport where transport occurs via diffusion of ligands bound to carrier molecules 
in the cell membrane 11, Q], see appendix iDl 

Our theoretical description of morphogen transport captures the key processes that are 
relevant for ligand transport and ligand-receptor interactions in multicellular epithelia such 
as the wing disk. We used a simplified description of these processes and neglected several 
aspects that could play a role. For example, we did not account for cell divisions and 



tissue growth in our description 



271 ]. More importantly, the presence of different receptor 



types which is quite common for signaling molecules of the TGF-/3 superfamily like Dpp 
was neglected in our description. These receptors typically form dimers or other complexes 
and, in general, the affinity of the ligand is different for the different receptor types and 
complexes. The trafficking of ligand-receptor pairs inside the cell is a very complex process 
that is only crudely captured in our description by a few parameters. 

Our coarse-graining procedure starts from a discrete cellular representation and allows 
us to obtain effective transport equations in a continuum limit. This provides a theoretical 
framework for a quantitative analysis of the spreading and trafficking of signaling molecules 
in and between cells. Using this approach one can relate experimental data obtained at 
different scales ranging from the cell to the tissue level. For example, the situation shown in 
Fig. [7| mimics recent experiments done in the Drosophila wing disk in which endocytosis is 
blocked in patches of cells in the tissue (see section [TTJ) [6] . The calculation results shown in 
ig.Elare consistent with the experimental data obtained in these shibire clone experiments 
They show a ligand depletion of decreasing relative magnitude ("contrast") behind 
the clone region as it is observed experimentally, see Fig. H. Note, that the clone region 
itself is devoid of ligands which is also seen experimentally and is evidence for a transport 
mechanism via transcytosis jlHI ]. 

Transport processes of signaling molecules in tissues show many features in common 
which are captured by our description. Our study highlights some general properties of 
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these systems such as the robustness of gradients which are largely independent of parameter 
values and molecular details. We expect that such general features can play an important 
role in very different biological signaling systems. 



APPENDIX A: DERIVATION OF THE EFFECTIVE TRANSPORT EQUA- 
TIONS IN ONE DIMENSION 

To perform a continuum limit of Eq. we introduce the densities l(t,x), r^^ T \t, x), 
ri(t, x), s^ r \t,x), and Si(t,x), such that x = na, L n (t)/a = l(t,x), R$(t)/a = 
r®(t,x), R®(t)/a = rM(t,x), E®/a = n{t,x), S^(t)/a = 8®(t,x), S^(t)/a = 
s^(t,x), and S$/a = Si(t,x). Kinetic equations for these are obtained by replacing the 
discrete densities R$, R^ r \ S$, S^ T \ and L n in Eq. (JTJ) with the continuum densities 
Tii r« r ', s ( - 1//r \ and I. The spatial separation of the quantities defined on the lattice 
as indicated in Fig. is taken into account by including terms up to second order in a 
power series expansion in x. We also change variables to r±(t,x) = r^ l \t,x) ± r^(t, x) and 
s±(t, x) = s^(t, x)±s^ r \t, x) so that r + and s + measure the total free and ligand bound sur- 
face receptor concentrations per cell and r_ and s_ the polarization of these concentrations 
on the cell surface respectively. This yields the continuum equations 

d t l = k oS s + - (ak on r + + e deg )Z - ab ^ on ld 2 x r+ + h J^l d 2 x s + + D d 2 x l 

o o 

abk m bk oS 

■ ld x r_ + — — d x s_ 



2 2 

n 1 1 , l n , ; \ , a b ^on o2r , \ a ~~ ®) ^ext 2 Clbk on 

d t s+ = ak on lr + + b cxt Si - (6 int + k oS )s + H — r + d x l H d x s { — r_d x l 



a - 6) 6mt o 2 (a - 6)6 



d t Si = ~(b cxt + b deg )si + b mt s + -\ d x s + d x s- 

2 

a li ru , i \ i a ^ ^ on a2i abk on (a — 6)6 ex t a 
d t s- = ak on lr_ - (b int + fc off )s_ H — r^dj — r + d x l H d x Si 

2 

d t r + = — + f ext ri - f- mt r + - ak on lr + + k oS s + - —r+d^l 
a 



2 . 

I ext r 



2 . 



d t n = -(/ cxt + f dcg ) ri + / int r + + (a g )/int d 2 x r + - (a 2 6)/int d x r. 

ab 2 k on 2 abk on (a - 6)/, 



d t r- = -f- mt r- - ak on lr„ + k oS s_ — r_d x l H — r + <9J H — - d x Ti. (Al) 

o Z Z 

Here, it is reasonable to neglect derivatives of higher order with respect to x because the 
most important contribution to ligand transport on large length scales comes from the second 
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derivative terms. This is due to the fact that first derivative terms must not appear in an 
effective transport equation because of the mirror symmetry of the original description. 

In the absence of degradation and production, there are two conserved quantities in the 
system, namely the total ligand number and the total receptor number. Indeed, the kinetic 
equation for the total ligand density A = I + Sj + s + that follows from Eq. (jAlj) can be 
written continuity equation with sink term 



d t X = -d x j - b dcg Si - e deg L (A2) 



Here, the total ligand current is 



abk on (a-b)b int bk oS ab 2 k on 
J = -^~ lr - + g S_ 2~ S ~ 8 — ~ + 

(a-b) 2 b cxt (a-b)\ nt + b 2 k oS 

-o x s t d x s+ - D d x l. (A3) 



8 8 
Note, that the terms involving r_ and s_ appear directly whereas all other terms are 
proportional to derivatives of r + , s + , s^, or I. The kinetics of the total receptor density 
p = V{ + r + + Si + s + is given by another continuity equation with source and sink terms: 

d t p = -d x j p + — - f deg ri - b deg Si (A4) 

(X 

with the total receptor current 

j P = a n b (fmtr- + b- mt s-) - ^ ^ (f cx td x ri + f- mt d x r+ + b ext d x Si + b int d x s + ) . (A5) 

Z o 
The individual terms of the currents j and j p are difficult to interpret. However, it will 
become clear below that the terms in j p do not give rise to transport over large distances 
whereas this is the case for the terms in j. 

The equations (|A2|) and (|A4|) have the unpleasant property that they relate the time 
development of A and p to that of all the individual quantities I, Sj, s±, Tj, r± whose time 
development is given by the set of coupled partial differential equations (|A1[) . It would be 
better if the kinetics of A and p could be described by equations which only involve these 
two quantities. This can be achieved by exploiting a separation of time scales. 

As discussed in the main text, the relaxation time scale r a for the kinetics in one cell is 
much smaller than the time scale t^ d for ligand transport on a large length scale £d on which 
the ligand profile develops. As we are interested in the behavior on large length scales in the 
continuum description, we exploit r a t^ d by making an adiabatic approximation in which 
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the system equilibrates infinitely fast locally. This is done by setting all time derivatives 
in Eq. (jAljl to zero and neglecting the second derivative terms. The resulting equations 
provide 5 relations between the 7 variables I, Tj, r±, Sj, s±. Here, we also assume that the 
production and degradation rates are small compared to the other rates. Formally, this 
corresponds to setting 6 deg = ed og = /dog = / sy n = 0. This procedure yields the relations ®. 

Note, that we have kept the first derivative terms for r_ and s_ in Eq. This is done 
to retain all second derivative terms when inserting Eq. (jSJ) into Eqs. (jA2j) . (jA3|) . (jA4|) . 
and (|A5|) . Using Eq. (jSJ), one can express /, r i; r_, s i; s_ in terms of r + and s + and spatial 
derivatives thereof. Finally, r + and s+ can be expressed in terms of p = r j + r + + Si + s + and 
A = / + Sj + s + . Mathematically, there exist two solutions for s+(A, p) and r + (A, p) but only 
one of them satisfies the physical requirement that s + (A = 0,p) = and r + (A,p = 0) = 0. 
Thus, we can uniquely express l,ri,r±,Si, and s± in terms of A and p in the adiabatic 
approximation. 

Using these expressions, it is straightforward to cast Eqs. (|A2|) . (|A3j) . (|A4jl . and (|A5Jl 
into the two coupled partial differential equations (jOJ) and 0. The explicit expressions for 
the transport and degradation coefficients in this one dimensional description of morphogen 
transport are 

— d 3 b ex tbi n t(b ex t + &int)/ext 2 /intfcoff^onP[ — 26i nt (6 ex t(/ext + /int)&off 

+a(b cxt + b- mt )f cxt k on X) 2 + 2v4(A,p)(-(6 int / cxt (2/ int (6 int + k oS ) + ab int k on X)) 
+b C xt(— 2/ ext / int A; off + &int(/int^off + /ext(~2/ int + fc D g — a& on A)))) 
+2a6 int (6 cx t + b int )f cxt k on (A(X,p) - 26 cxt (/ cxt + f- m t)k oS 
+2a(6 ex t + bi Qt )f cxt k on X)p - 2a 2 6 in t(6 C xt + ^nOVerf^anV] -1 

+2Z) a6 ext (6 ext + 6int)/ext(/ C xt + /int)^off^onP[^(A,p)(v4(A,p) 
-&ext(/intfcoff + /extOoff + a^on(A - p))) + ab mt f cxt k on (p - A))]" 1 

D p (X,p) = a% xt b in t(b ext + &int)/ext 2 /intfcofffcon 2 Ap[2((-(/i nt (6 int + fc off )) - a& int £; on A) 
(&ext(/ext + fmt)k $ + a(b ext + b int ) f ext k OI1 \) (A(X, p) 
+b cxt (fcxt + fintjkoB + a(b cxt + bmt) f cxt k on X) 
+a(b cxt + b- mt )f eXk k on (-(A(X,p)f int (b int + k oS )) + aA(X, p)b int k on X 

+2(/in t (&int + &off) + a&int^onA)(-(6 xt(/ext + /int)&off) 

-\-a(b cxt + b[ nt ) f ext k OQ X)) p — a 2 (6 cx t + feint) fext 2 k on 2 

x(/int(&int + ^off) + a&intfconA)p 2 )] _1 
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+2D b cxt (f cxt + /int)fcoff((((/e3rt + frntjhs + af cxt k on X) 2 

+«/ext^on((/cxt + /int)&off — o/ext^onA)p)fe e xt 2 

-A(A,p)((/ cxt + f int )k oS + af cxt k on X)b cxt + ab mt f cxt k on (2f cxt X{k oE + a£; on A) 
+/ C xt(^off - 2ak on X)p + f int k oS (2X + p))b cxt + ab int f cxt k on X 
x(a6 int /cxtA;on(A - p) - A(A,p)))[A(A,p)(,4(A,p) - b cxt (f int k oS 
+fe X t(k oS + a£; on (A - p))) + ab mt f cxt k on (p - A)) 2 ] -1 

fc A (A, p) = [&d eg &int(a&int/cxt/i;on(A + p) " A(A, p) + b cxt (f int k oS 

+ /ext(&off + afc on (A + p))))][2a(6 ext + feintf/ex^onA]" 1 
fcp(A, p) = [(- ((feeaet + &int) /deg/int) + & deg "int (/ex t + /int))(-A(A,p) 

+^cxt (/eart + /int) ^off + O (&ext + &int) /ext^onA) 

(&ext + &int) /eart ((&ext + &int) /deg/int + ^deg&int (/eart + /int)) &onP] 

X [2a(6 cxt + 6 int ) 2 /ext (/ext + /int) ^onP] ~* , with 



A(X, p) = [(6 ex t(/int^off + fext(k s + ak on (X - p))) + ab int f ext k on (X - p)) 2 

+4a6 ex t(&cxt + &int) /ext (/ext + Znt)&off&onP] 1/2 - (A6) 



APPENDIX B: DIRECTIONAL BIAS OF INTRACELLULAR TRAFFICKING 

We study the effects of a bias in the description of transcytosis. We introduce a di- 
mensionless parameter —1/2 < (3 < 1/2 which measures this bias in Eq. ([!}. The kinetic 
equations of the discrete description with bias are: 

j t S { ; ] = -k oS S^ + k on R^L n -b int S^ + b cxt (l/2 + P)S® 

j t S® = -k oS S® + k on R®L n ^ - btoS® + &ext(l/2 - (3) S®, (Bl) 

with the kinetics for the remaining quantities as in Eq. (0). For (3 = —1/2, all receptor- 
bound ligands are externalized on the left surface, for (3 = 1/2 on the right surface, and for 
(3 = 0, we recover the unbiased description. This implies that the externalization of the free 
receptors remains unbiased. We proceed as before to derive the continuum equations. The 
adiabatic approximation changes to 



k on alr_ — 2(3b ex \ i Si — ab! ^ on r + d x l + - — ^ 9 x s 



(B2) 
bmt + k oS 
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with the other relations as in Eq. (|5|). Finally, the transport equations (JHJ) and (JJJ) generalize 
to Eq. Q with the new effective drift velocity 

Vp{\,p) = 2a 2 pb cxt b int f cxt fi nt k oS k on p[b cxt (f illt (b int + k oS ) + ab int k OQ X)((f cxt + f int )k s 
+af ex tk on X) + ab ext f ex tkon(fmt(bmt + k oS ) - ab int k on X)p 
+(/int (fat + ^ofr) + ab- mt k on X)A(X, p) + ab- mt f cxt k on (f- mt k o{i (\ + p) 
+b- mt {ak on X{X -p) + f int (X + p)))]-\ (B3) 

where A(X, p) is defined in Eq. (IA6|) . The other coefficients appearing in Eq. © remain as 
in the case without directional bias. 



APPENDIX C: SIMPLE LIMITS OF THE EFFECTIVE TRANSPORT EQUA- 
TIONS 

In biological systems, it is well possible that some of the processes included in our de- 
scription of ligand transport are much faster than the others. For example, the binding of 
the ligand to its receptor can be fast compared to other processes due to the small volume 
of the gaps between cells. The confinement of the ligands to this small volume leads to 
frequent collisions between ligands and receptors. Assuming that the reaction is diffusion 
limited this can lead to a high reaction rate k on . 

It is worthwhile to note that Eqs. © and (|7J) become simpler if this or another one 
of the cellular processes is much faster than the others, i.e. if the corresponding rates in 
our description are very large. For example, in the limit of very fast binding of ligands to 
receptors, i.e. for k OD — > oo, we find 

^ / , \ ~ a2 b C xtbintfextfmtkoSP 

Da(a ' p) = b(M 

j — x / \ \ o b cx tb- mt f ex tf m tk fi\ 
D p {X,p) = — 



B(X, P ) 

h, h 

k\{\p) 



bdcg feint 



where 



feext + feint 

, /, \ fedegfeint(/cxt + /int)A + (fecxt + feint ) /dcg/int (ft ~ X) 

A ' Pj " (feext + feint) (/ext + /int)p ' 1 ' 



B\\P) — 4(feint/ext/int (feint + ^off) + fecxt (feint/cxt/int + /ext/int^off 

-femt(/ext + fmt)k oS ))X - 4(6 ext + fe int ) /ext /int (feint + k oS )p. 
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There are no free ligands in this limit because I = via Eq. (J3J). As all ligands are bound to 
receptors, free diffusion does not contribute to the current and Dq does not appear in Eq. 
(|C1|) . This also imposes the constraint A < p. The example given by Eq. (|C1|) is instructive 
because the effective transport and degradation coefficients are much simpler than those in 
the general case (jA6|) . Many of the properties discussed above for the general case that is 
valid for arbitrary k on can be read directly from the expressions in Eq. (|C1|) . For example, 
if b ext = ligands do not move because D\ = D p = in Eq. (jCf|) . Due to the constraint 
A < p, however, the statements for the asymptotic behavior for A — > oo do not apply to 
this case. Furthermore, D\ does no longer exhibit a maximum as a function of A. It either 
grows or decreases monotonically depending on the parameter choice. 

In principle, one can write down simpler expressions as in Eq. (jCIj) for many different 
limits. If several transport steps are much faster than the others, only the ratios of the 
corresponding parameters enter the simplified description. For example, if the binding and 
un-binding of the ligand to the receptor is faster than all other processes the effective diffusion 
coefficient and degradation rate do not depend on k on and k s individually but only on 
the ratio k on /k ^. The number of parameters can thus be reduced to obtain the minimal 
description for a given situation. A simple but instructive example is the situation k on ^> 
Ks 3> b cxt 3> &int in Eq. (fTTj). for which we obtain the effective diffusion coefficient 

lim lim lim -D(A) = 6; nt a 2 /4. 

frcxt^OO fcrff-^OO fcon^OO 

This reflects that the only slow process — in this case the internalization of ligands at rate 
6 int — limits the transport efficiency and defines the effective diffusion coefficient. Note, 
that the limits taken above do not commute because D(\) = for k g ^> k on . 

APPENDIX D: TRANSPORT BY DIFFUSION IN THE CELL MEMBRANE 

As another example for an application of our theoretical framework, consider a transport 
mechanism in which ligands can move across cells by passive diffusion in the cell membrane. 
The ligand first binds to a receptor molecule on the cell surface. This complex then diffuses 
in the cell membrane. At any time, the ligand can detach from the receptor it occupies 
and, after diffusing over a short distance in the extracellular space, it can attach to a new 
receptor that can be located on the surface of the same or a different cell. Note, that if we 
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simply replace the receptors with HSPG molecules this trans por t mechanism is very similar 



111- 



to one that was recently suggested for Dpp in the wing disk 

We can describe this mechanism on the same lattice structure as used above for transcy- 
tosis, see Fig. ^ For simplicity, we focus on the essence of the transport phenomenon and 
do not include ligand degradation or the production and degradation of receptors in our 
description. Furthermore, the receptor concentration is assumed to be constant on the cell 
surface. This description is similar to that of transcytosis with constant surface receptor 
concentration given by Eq. (|10p. However, instead of the internalization and externalization 
of receptor-ligand complexes as in transcytosis, dislocation of the ligand-receptor complexes 
across one cell is due to diffusion on the cell surface. In a discrete one dimensional descrip- 
tion, this effect is captured by hopping at a rate h ~ D M /2a 2 between the left and right 
surface of each cell, where Dm is the diffusion coefficient of receptors in the cell membrane. 
Using the same notation as in Eq. (fTUj) . a discrete description of this mechanism reads 

— L n = k s(S$ + S n+ i) — k on (R — — Sn + i)L n 

= _ kosS M + koa ( R/2 _ S W)L n + h(S® - S£>) 
j t S^ = _ ko{fS U + k on {R/2 - 5«)L n _! + h(S® - 

Applying the same method as in section IIIIBl we obtain the continuum transport equa- 
tion dfX = d x (Dsr>Wd x X) where A is the total ligand concentration and 

a 3 rhk oS k on (E(X) - k oS + ak on (-A + r)) 
SDl j ~ AE(X)(E(X)h-h(k oS + ak on X) + ar(h + k oS )k on ) ,W1 

E(X) = [-4a 2 rA; on 2 A + (A; off + aA; on (r + A)) 2 ] 1/2 . (Dl) 

Dsd(X) has properties that are similar to those of the effective diffusion coefficient of the 
transcytosis model. It exhibits a maximum as a function of A and decays as Z^sd ~ A -2 for 
large A. Hence, upon inclusion of ligand degradation in the description, gradient formation 
by this mechanism exhibits similar properties as transcytosis. 



APPENDIX E: TRANSPORT BY EXTRACELLULAR DIFFUSION 

In order to ultimately identify the mechanism of morphogen gradient formation present 
in a given system, it is important to develop mathematical descriptions of all potentially 
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relevant transport mechanisms so that these can be compared to the available experimental 
data. Extracellular diffusion is widely believed to be the dominant transport mechanism 
for some morphogens si, Q]. For this reason, we briefly discuss the extracellular diffusion 
dominated limit of the transport phenomenon defined by Eq. (JTJ in this section. 

Assuming that extracellular diffusion gives the dominant contribution to the ligand cur- 
rent, we neglect the contribution of transcytosis. This approach is valid if the rates of 
ligand trafficking are such that the effective diffusion coefficient resulting from transcytosis 
D(X,p) \ p= o in Eq. © is much smaller for all values of A and p than the extracellular diffusion 
constant D . 

The derivation of the continuum limit of Eq. (JTJ) is straightforward in this situation. 
Only the equation for the time development of the free ligand concentration contains a 
linear diffusion term. All other quantities can be described by a coupled set of ordinary 
differential equations. Together, these constitute a set of react ion- diffusion equations. With 
the notation used throughout this article, the kinetic equations in two dimensions read 

d t l 
d t s+ 
d t Si 
d t r + 



D Al - k on a 2 lr + + k oS s + - 
k on a 2 lr + - (b mt + k oS ) s + + & cxt S; 

&int s + — {bext + ^deg) s i 

^ + k oS s+ + /extn - k on lr + - f int r + 



/mt?"+ — (fext + /deg) r i- 



This is essentially the system that was previously studied in 



(El) 
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Jjj . This reaction-diffusion 



system can reproduce some experimental results that were obtained for the morphogen Dpp 



and its receptor Tkv in the wing disk u&Q. However, there are experimental observa- 
tions that are in disagreement with the solutions of Eq. (|E1|) indicatingthat this is not 



a correct description of the ligand and receptor kinetics in this system 
experimental systen 
during development 



131 ] . In other 



experimental systems, reaction-diffusion mechanisms may play a role in pattern formation 
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FIG. 1: Schematic representation of ligand transport by transcytosis in a chain of cells of diameter 
a indexed by n. The rates of ligand-receptor binding and un-binding, internalization and external- 
ization of ligand-receptor pairs are denoted k on , k Q s, 6; n t and 6 ext . Degradation of ligand occurs 
inside the cells with rate &deg and in the extracellular space with rate edeg- Ligands can also hop 
directly between neighboring extracellular spaces at a rate 2Do/a 2 which describes their movement 
in the extracellular space around the cells by passive diffusion with diffusion coefficient Dq (not 



shown). Figure modified from 



2C|. 



FIG. 2: Time development of gradient formation in our description of ligand transport. Ligand 
densities in the presence of a source at x = at different times tb deg = 0.72, 2.16, 3.6 during gradient 
formation (black lines) and in steady state (red lines). Lines indicate solutions to Eq. (|6"J). while 
symbols indicate solutions to Eq. (JJ) for comparison. (A-C) Time development of the profiles of the 
total ligand density X(x, t) (A), the total receptor density p(x, t) (B), and the receptor bound ligand 
density Si(x,t) + s + (x,t) (C) in the absence of extracellular diffusion, i.e. for Dq = 0. (D) Like 
A but with 2D Q /a 2 b deg = 10/3, i.e. in the presence of extracellular diffusion. All concentrations 
are normalized to the steady state value of the surface receptor concentration in the absence of 
ligands tq. Initial conditions at t = 0: A(x) = and p(x) = (1 + /int//ext)^o- Parameters are: 
Ks/bdcg = b int /b dcg = /int/^dcg = 1000/3, k on ar max /b deg = 8000/3, b cxt /b dcg = f e *t/b deg = 2000/3, 
edeg/frdeg = 2/3, f deg /b deg = 1, / s ° yn /ar max 6 deg = 1/12, ip = 2, jo/b deg ar mSliX = 25/6, r /r max = 1/7, 
and j = at x/a = 50. 
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FIG. 3: Effective transport coefficients and degradation rates for transcytosis. (A,B) Coefficients 
D\(\,p), D p (X,p), and k\(X, p) in the transport equation © as a function of the dimensionless 
ratio X/p of the total ligand concentration A and the total receptor concentration p. The solid lines 
show the coefficients D\ and D p in presence of extracellular diffusion with 2Do/° 2 ^deg = 10/3. 
Broken lines show these coefficients in the absence of extracellular diffusion, i.e. with Dq = 0. 
Inset in A: the solid line shows k\ with edeg/^deg = 2/3 and the broken line with ed C g = 0. 
po denotes the steady state total receptor concentration in the absence of ligands. (C) Effective 
diffusion coefficient D(X) in the transport equation for the constant surface receptor approximation 
as a function of the ligand concentration X/r for 2Do/° 2 &deg = 10/3 (solid line) and Dq = 
(dashed line). Inset: effective degradation rate k(X) as a function of X/r for ed eg = (dashed line) 
and edeg/&deg = 2/3 (solid line). The constant total surface receptor concentration is denoted r. 
Parameters as in Fig. [2] with k on ap/bd e g = 8000/3 in A,B and fc n£w/frdeg = 8000/3 in C. 



FIG. 4: Ligand transport by transcytosis with a directional bias. (A) Drift velocity Vg from 
Eq. ()B3|) as a function of X/p. (B,C) Time development of gradient formation with directional 
bias. Profiles of the total ligand concentration X(x, t) (B) and the total receptor concentration 
p(x,t) (C) in the presence of a source at x = at different times tftdeg = 0.72,2.16,3.6 during 
gradient formation (black lines) and in steady state (red lines) . Lines indicate solutions to Eq. |2| , 
while symbols indicate solutions to Eq. (|B1|) for comparison. All concentrations are normalized 
to the steady state value of the surface receptor concentration in the absence of ligands ro- Initial 
condition: X(x) = and p(x) = (1 + /int//ext)^o- Parameters as in Fig. Elwith (3 = .1, Dq = 0, and 
j = at x/a = 100. 
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FIG. 5: Non-robust and robust steady state gradients in our description of morphogen transport 
with constant surface receptor concentration given by Eq. (|llj) . (A) Ligand profiles in the steady 
state for Dq = and jo/b^R = 7 where robustness is small, 1Z ~ 0.1. The profile is strongly 
affected by halving (dotted line) or doubling (dashed line) the ligand current of the reference 
state (solid line). The positions of an arbitrarily chosen concentration threshold are indicated. 
(B) Ligand profiles in the steady state for Dq = and jo/b deg R = 70 where robustness is large, 
1Z ~470. The reference profile (solid line) is almost unaffected by halving (dotted line) or doubling 
(dashed line, covered by the solid line) the ligand current of the reference state. (C) Like B but 
with extracellular diffusion (Do/a 2 bd eg = 50) which reduces robustness to 71 ~ 0.32. The insets 
in A-C show the respective profiles of the receptor-bound ligand concentration s+ + Sj which is 
the biologically relevant quantity that triggers signal transduction in the cells. (D) Robustness 1Z 
of steady state ligand profiles as a function of the ligand current j$ from the source for different 



values of the ratio of the extracellular diffusion length ^ = y-Do/ e deg and the cell size a. The 
description of morphogen transport with constant surface receptor approximation given by Eq. Qll|l 




was used to calculate 1Z. Figure modified from 20]. Parameters are &i n t/&deg 




3x 10 3 



k on ar/b de g = 1.1 X 10 4 , ed cg /&dcg = 5, and k e/b de g = 7 x 10 2 . 
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FIG. 6: Ligand transport in two dimensions. (A) Tissue in the region of the wing disk where the 
Dpp gradient forms. Cell membranes are labeled in red, the morphogen Dpp is shown in green. (B) 
Schematic of this tissue. Source cells which produce ligands are shown in darker gray, the cells in 
the target tissue in light gray. (C) The rates for the various cellular processes are denoted as in one 
dimension, see Fig. ^ For the triangular lattice with hexagonal cells which we use here, receptor 
bound ligands can be present on the six edges of each cell and inside the cells. The concentration 
of the receptor-bound ligands on edge j of cell n is denoted S n j, that inside cell n is termed S$ ■ 
Free ligands exist in the gaps between two cell surfaces. Their concentration in gap (n,j) which 
is located adjacent to edge j of cell n is denoted L/ n j\. (D) Triangular lattice structure with 
hexagonal cells used in our discrete theoretical description of ligand transport in two dimensions. 
This lattice structure approximates the situation shown in B. The source cells which are shown in 
darker gray secrete ligands into the extracellular spaces surrounding them. The cell diameter is a. 
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FIG. 7: Time development of ligand densities X(x, t) in two dimensions. The solution to the 
continuum transport equation (I22|) is compared to that of the discrete description given by Eq. 1)21(1 
in presence of a region which the ligand cannot enter located at 6 < x/a < 11 and —4 < y/a < 4. 
In the discrete description, this is realized by setting 6; n t to zero in the region. In the continuum 
description, zero flux boundary conditions are imposed on the outlining edge of the region. (A-D) 
Two dimensional ligand profiles at tfrdeg = 3.6 (A,B) and at t&deg = 17.3 which is close to the 
steady state (C,D). These profiles were obtained by solving the discrete (A,C) and continuum 
description (B,D) respectively. The rectangular region which the ligand cannot enter appears in 
white. It appears smaller in A and C than in B and D because in the discrete description, ligands 
can still bind to surface receptors on the cells located at the edge of the region so that A > for 
these cells. (E-G) Profiles of \(x,t) along the slices indicated in B at i^deg = 0.72,2.16,3.6 (black 
lines) and in steady state (red lines). Lines indicate solutions to Eq. (j22j), while symbols indicate 
solutions to Eq. (|21|) for comparison. (H) Contrast of the depletion of \(x,t) shown in G. The 
contrast is defined as c(t) = 1 — Xb(t) / X (t) where \b(t) and A G (t) are the total ligand concentration 
at the locations shown by the crosses in A, i.e. directly behind the clone and far away from it, 
respectively. Initial condition: X(x) = 0. Parameters as in Fig. [2]with Dq = 0, io/&dogO r = 25/3 
at x = 0, j = at y/a = ±25 and at x/a = 50. 
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